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ABSTRACT 

We report on the discovery of 10 additional galaxy clusters detected in the 
ongoing Swift/BAT all-sky survey. Among the newly BAT-discovered clusters 
there are: Bullet, Abell 85, Norma, and PKS 0745-19. Norma is the only cluster, 
among those presented here, which is resolved by BAT. For all the clusters we 
perform a detailed spectral analysis using XMM-Newton and Swift/BAT data 
to investigate the presence of a hard (non-thermal) X-ray excess. We find that 
in most cases the clusters' emission in the 0.3-200 keV band can be explained 
by a multi-temperature thermal model confirming our previous results. For two 
clusters (Bullet and Abell 3667) we find evidence for the presence of a hard X-ray 
excess. In the case of the Bullet cluster, our analysis confirms the presence of a 
non-thermal, power-law like, component with a 20-100 keV flux of 3.4x10"^^ erg 
cm^^ s^^ as detected in previous studies. For Abell 3667 the excess emission can 
be successfully modeled as a hot component (kT=~13keV). We thus conclude 
that the hard X-ray emission from galaxy clusters (except the Bullet) has most 
likely thermal origin. 
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Subject headings: galaxies: clusters: general - acceleration of particles - radiation 
mechanisms: non-thermal - magnetic fields - X-rays: general 

1. Introduction 

The study of clusters of galaxies at X-ray energies is key to understand the mechanisms 
that heat the intra-cluster medium (ICM) and to measure the pressure due to cosmic rays 
(CRs), magnetic fields and turbulence. In particular shock heating can be infiuenced by CRs 
if a significant part of the shock energy is transferred to charged particles. Indeed, large-scale 
shocks that form during the process of cluster formation are believed to be efficient particle 
accelerators (e.g. Sarazin 1999; Ryu & Kang 2003). Thus the pressure support of CRs to 
the ICM might be relevant. 

Not surprisingly, the role of CRs in the formation and evolution of clusters of galaxies 
has been much debated. Churazov et al. (2008) suggest that in massive galaxy clusters 
hydrostatic equilibrium is satisfied reasonably well, as long as the source has not experienced 
a recent major merger. On the other hand, other studies (e.g. Miralda-Escude & Babul 1995; 
Nagai et al. 2007) showed that the non-thermal pressure due to CRs, magnetic fields and 
micro-turbulence can affect the mass estimates based on hydrostatic equilibrium. Knowing 
the importance of CRs, the mechanisms that heat the ICM and the frequency at which it is 
shocked, is crucial for the upcoming X-ray and Sunyaev-Zeldovich effect cluster surveys (see 
Ando & Nagai 2008). 

The detections of extended synchrotron radio emissions (e.g. Willson 1970; Harris & 
Miley 1978; Giovannini et al. 1993; Giovannini & Feretti 2000; Kempner & Sarazin 2001; 
Thierbach et al. 2003) represent the main evidence that a population of non-thermal rela- 

tivistic electrons exists in the ICM. These very same electrons can produce X-rays via Inverse 
Compton (IC) scattering off cosmic microwave background (CMB) photons (e.g. Rephaeli 
1979; Sarazin 1999), or via non-thermal bremsstrahlung (e.g. Sarazin 1999; Sarazin & Kemp- 
ner 2000) or synchrotron radiation (Timokhin et al. 2004; Inoue et al. 2005). Detecting this 
non-thermal radiation is difficult because of the bright and dominant ICM thermal emission. 
Studying clusters above 15keV, where the intensity of the thermal component decreases 
quickly, might prove to be an effective probe of the non-thermal emission processes. In- 
deed, in the past, the detection of non-thermal emission in the hard X-ray spectra of a 
few galaxy clusters has been reported (see e.g. Kaastra et al. 2008, for a complete review). 
However its actual presence and origin remain controversial (e.g. Rephaeli et al. 1987, 1999; 
Rephaefi & Gruber 2002; Rossetti & Molendi 2004; Sanders et al. 2005; Renaud et al. 2006; 
Fusco-Femiano et al. 2007; Lutovinov et al. 2008; Molendi & Gastaldello 2008). 
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In a first paper (Ajello et al. 2009a), we reported about the detailed analysis of 10 
galaxy clusters serendipitously detected in the Swift/BAT all-sky survey above 15keV. In 
that study we concluded that there were no significant evidences for the existence of hard 
X-ray excesses detected in the spectra of clusters above the BAT sensitivity. In this paper 
we report the analysis of 10 additional clusters that have been recently detected, thanks 
to the deeper exposure, in the ongoing BAT survey. We combine BAT and XMM-Newton 
data to find the best spectral fit. Assuming that there is a non-thermal emission due to IC 
scattering on CMB photons, we estimate the upper limit of its fiux in the 50-100 kev band. 
This information allows us to estimate the intensity of the magnetic fields in these galaxy 
clusters. 

Throughout this paper we adopt a Hubble constant of Hq — 70kms~^ Mpc~^, Q.m = 
0.3 and Qa — 0.7. Unless otherwise stated errors are quoted at the 90% confidence level 
(CL) for one interesting parameter and solar abundances are determined using the meteoritic 
values provided in Anders & Grevesse (1989). 

2. Clusters in the Swift/B AT Survey 

The Burst Alert Telescope (BAT; Barthelmy et al. 2005), on board the Swift satellite 
(Gehrels et al. 2004), represents a major improvement in sensitivity for imaging of the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° 
partially coded) sensitive in the 15-200 keV domain. As shown in several works (e.g. Ajello 
et al. 2008a, 2009b; Tueller et al. 2010), thanks to the deep exposure, BAT reaches sub- 
mCrab (e.g. < 10^^^ erg cm~^ s~^) sensitivities in the entire high-latitude sky. This already 
allowed BAT to detect 10 galaxy clusters above 15kcV after ~2 years of all-sky exposure (see 
Ajello et al. 2009a, for details). With ~6 years of all-sky exposure acquired, the sensitivity 
of BAT has increased substantially leading to the detection of almost 1000 source in the 
hard X-ray sky (see Cusumano et al. 2010, for details). Here we present a detailed spectral 
analysis of 10 galaxy clusters detected in the BAT survey of Cusumano et al. (2010). 

2.1. Analysis of Swift/BAT Data 

The BAT makes images of the sky thanks to a coded mask (with a random pattern) 
placed above a position sensitive detector plane (see Barthelmy et al. 2005, for details). The 
sky radiation passing through the aperture is coded by the mask pattern and recorded in the 
detector plane. The pattern of the mask is such that a source at a given position in the FOV 
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casts a unique shadow onto the detector plane and thus its emission can be easily deconvolved. 
The randomness of the mask pattern ensures that the cross-talk between sources (e.g. some 
flux from a given source is wrongly attributed to another one) at different positions in the 
FOV is minimum. Moreover, to minimize this and other systematic uncertainties that can 
arise in the BAT survey, Swift adopts a random roll-angle strategy when pointing at the 
same position in the sky. This means that whenever Swift is pointing at a given direction 
in the sky, the roll angle (e.g. the angle on the plane orthogonal to the pointing direction) 
is chosen randomly within a range of the nominal (e.g. within ±2 degrees) pointing. This 
ensures that pointings are never exactly the same and that sources never fall in the same 
relative positions in the BAT FOV. 

A decoding procedure is required in order to reconstruct the original sky image. A 
variety of methods can be used to reconstruct the sky image in the case of a coded mask 
aperture (see Skinner et al. 1987; Ajello et al. 2008a, for a general discussion on reconstruction 
methods). Among them, standard cross-correlation of the shadowgram (e.g. the information 
recorded on the detector plane) with a deconvolution array, the mask pattern, via fast 
Fourier transforms (FFTs), is the most often used. Generally, sky images are obtained 
for each individual observation, where an observation is defined as a period during which 
the attitude is stable and constant. Subsequently, another procedure, such as resampling 
and re-projecting, is needed in order to assemble the final all-sky image. As discussed in 
the Appendix in Ajello et al. (2008c), the balanced correlation (Fenimore & Cannon 1978) 
used to deconvolve BAT observations and source spectra performs a standard background 
subtraction (see the above references for the exact implementation of this method). However, 
this technique works well as long as the background in the array is fiat and not correlated 
with the mask pattern. The background in the BAT array is not fiat due to the presence 
of many background components, the brightest of which is (below '^60keV) the Cosmic 
X-ray Background (sec Ajello ct al. 2008b, for details). Thus, the balanced correlation 
alone provides imperfect results and produces a noticeable background contamination in 
the sky observations and in the source spectra. This background contamination has been 
estimated (see Ajello et al. 2008c) to be <2 % of the Crab Nebula intensity in the 14-195 keV 
band (e.g. the BAT band). Thus, this contamination does not pose problems for strong 
sources, but becomes very relevant for the (spectral) analysis of faint objects with ~mCrab 
intensities as the clusters of this work. In order to correct for this residual background 
contamination we use the recipe presented in Ajello et al. (2008c). We use several templates^ 
of the BAT background (for each channel) which are fit together with the contribution of 



^For reference see the description of the batclean tool available at 
http://heasarc.nasa.gov/lheasoft/ftools/headas/batclean.html. 
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point sources to the BAT detector counts. In each observation the residuals are analyzed to 
check for additional structures and deviation from Gaussian statistics. If those are found, 
then thousands of residual maps are averaged together (in image coordinates) to create a 
blank field ohservaiion^ . These blank field observations become part of our template library 
of background models and fit once again to any observations which are being used. The 
process of residuals-inspection and template-creation is repeated until the residuals do not 
show any systematic feature. For a given observation, the last template model to be added 
is generated from observations which are close to it in time. This ensures that long-term 
variation in the BAT background (due e.g. to the orbit, activation of the spacecraft, noisy 
pixels etc.) are correctly taken into account. 

Adopting this technique and filtering the data in the way described in details in Ajello 
et al. (2008a), we extracted a 15-195 keV spectrum for the 10 galaxy clusters of this analysis. 
We used all the available observations at the time of this work (approximatively from 2005 to 
March 2010) resulting in an average exposure at each of the 10 positions larger than 14 Ms. 
It has already been shown that spectra extracted with this method are reliable and accurate 
over the entire energy range (see § 2.3 of Ajello et al. 2009a, for details). 

2.2. Studying "Extended" Sources with Swzft/BAT 

Coded mask telescopes are designed and optimized for the study of point-like sources 
(e.g. see Ajello et al. 2008b, and references therein). Formally the mask acts as a filter, 
canceling out those signals (celestial and not) whose spatial frequency is larger than the 
spatial frequencies of the mask tiles. This means that if an X-ray source extends over an 
area which is larger than the projection of the mask tile on the sky (i.e., 22.4' or, which is 
the same, the Full- Width at Half Maximum, FWHM, of the BAT Point Spread Function, 
PSF) then part of the X-ray flux is necessarily lost in the background. In simpler words, if a 
source is extended, then its shadow (produced when the radiation passes through the mask) 
on the detector plane looses part of its contrast (i.e., the mask is illuminated from all sides). 
The hmiting case is represented by the Cosmic X-ray Background, which extending over the 
entire sky, is completely removed by the BAT mask. Indeed in order to measure the Cosmic 
Background, different, non-standard, techniques have to be used (Ajello et al. 2008b). 

Clusters of galaxies are X-ray sources extending up to ~1 degree. In our first work Ajello 
et al. (2009a), we showed that all the clusters detected by BAT are consistent with being 



^Working in detector coordinates ensures that the contribution of any unsubtracted point source is aver- 
aged out. 
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point-like sources (for BAT) with the exception of Coma which is clearly resolved. Given 
the properties of coded masks expressed above, it is a good thing that clusters are seen as 
point-like sources. Indeed, in this case the flux measurement is correct while it is not if the 
source is detected as extended. 

In order to quantify this effect in more detail we performed several Monte Carlo simula- 
tions. We simulated an extended source whose surface brightness profile can be approximated 
by a beta model of the form: 



where is the core radius and /3 is typically in the 0.5-0.8 range (see e.g. Ettori et al. 
1998; Ajello et al. 2009a). We then reconstructed the sky image and detected the simulated 
source measuring its flux and the PSF FWHM. Simulations were performed, in the 15-55 keV 
band^, for several core radii and for a fixed ^ of 0.6. The results of this exercise are reported 
in Fig. 1. It is clear that a substantial part of the source fiux is lost if clusters have a core 
radius larger than 5' . This effect is a function of the emissivity profile of the source and 
depends on it. If we were considering a source with a Gaussian emission profile or a uniform 
profile (e.g. a disk-like emission) than the flux suppression would be even stronger. The 
presence of a cool-core acts in the opposite direction, indeed in this case a large part of 
the cluster's emission is confined within the inner few arc-minutes from the core (normally 
within the core radius) where the flux suppression (due to source extension) is negligible. 
Moreover, Fig. 1 gives a powerful tool to understand when the flux suppression takes place. 
Indeed, in all these cases the FWHM of the source is larger than the point-hke one (i.e. 
22.4' ). Coma which is resolved by BAT can be used as an example. Indeed, adopting a 
core radius of Tc = 10.7' (see e.g. Lutovinov et al. 2008), we can estimate from Fig. 1 that 
the measured FWHM should be ~ 26' . This is found in good agreement with the results 
reported in Ajello et al. (2009a). Moreover from the same graph we can estimate that the 
flux suppression is ~25 %. The next most extended clusters reported in Ajello et al. (2009a) 
are Ophiucus and Perseus. Ophiucus with a core radius of 3.2' (Watanabe et al. 2001) is not 
resolved by BAT and thus the flux suppression is negligible. Perseus has a core radius of 4.7' 
(Ettori et al. 1998) and thus at the limit where the fiux suppression might start playing a 
role. However, Perseus has a bright cool core whose emission profile can be modeled with a 
power law (e.g. see Ettori et al. 1998, and references therein). This concentrate most of the 
cluster's emission in the core which is never resolved by BAT. Hence, Perseus is detected as a 
point-hke object. The majority of the clusters detected by BAT, either those reported here or 



^The results of this analysis do not change if a different band is chosen. 
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those described in Ajello et al. (2009a), are detected as point- like objects. Norma (reported 
in this work) with a core radius of ~10' hke Coma, is also resolved by BAT. However, its 
marked North- West elongation and the presence of a bright nearby AGN make its case more 
complex that the simple spherically symmetric case discussed here. This, will discussed in 
details in § 3.8. 

To summarize this section, wc note that coded-mask telescopes are optimized for the 
study of point-like objects and can be used for the study of (intrinsically) extended objects 
only if these are detected as point-like (e.g. for a beta model this is true if the core radius is 
less than 5' ). If an object is detected as extended, then part of its flux (depending on the 
source brightness profile) is suppressed. The deviation of the PSF from the point-like PSF 
can be used to understand whether this effect is present. Among all the clusters detected 
by BAT, Coma and Norma (see next sections) are the only two sources which are extended 
in BAT and part of the source flux is lost in the BAT background. All the other clusters 
detected by BAT are not resolved by BAT and BAT can be safely used for the study of their 
emission. 



2.3. Analysis of XMM-Newton Data 

For all the clusters, we extracted a 0.5-8.0 keV spectrum using publicly available XMM- 
Newton (EPIC-PN) observations. The details of these observations are reported in Tab. 1. 

For each observation, XMM-Newton data are screened filtering for periods of flaring 
background. This is done examining the lightcurve in the 10-12 keV band and determining 
the rate of the quiescent background (the rate outside of the flaring episodes). For all the 
observations reported in this analysis this was very close to 0.5 counts s~^ in agreement with 
e.g. Nevalaincn et al. (2005) (and references therein). However, this is not sufficient for 
filtering out fiaring episodes that produce a soft background component (e.g. Nevalainen 
et al. 2005; Carter & Read 2007). We thus inspect the lightcurve in the 1-5 kcV band 
extracted in an annulus of inner radius 10' and outer radius 12' . We filter out all those 
times bins that deviate more than 3cr from the average quiescent background (determined 
through a Gaussian fit to the histograms of the rates) . 

We extracted the cluster spectrum using a single extraction region with a radius of ~10' 
. This is partly motivated by: i) the extent of the XMM-Newton EPlC-PN CCD (~12' in 
radius), ii) the fact that we use the two outer arcminutes to perform background filtering 
and in) the fact that for a typical beta profile (e.g. core radius of 3.8' and /3=0.7) this 
selection includes up to 94 % of the entire cluster emission and more than that if the cluster 
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Fig. 1. — Suppression of the source flux (solid line) due to the source extension as a function 
of core radius (rc) for an emission profile that follows a beta model (see § 2.2 for details). 
Error bars are of statistical origin and reflects the uncertainty in the reconstructed quantities. 
The dashed line shows the fractional increase of the FWHM of the source PSF (with respect 
to the PSF for a point-like source) when the source is "resolved" by BAT. For a cluster with 
a core radius of 30' the FWHM is ~30 % larger than the point-like source FWHM. 



has a cool core. The fact that the BAT PSF is consistent with the point-like one for all 
the clusters in this analysis (except Norma) suggests that the hard X-ray emission (above 
15keV) is coming from the inner part of the cluster. If this were not the case, we would have 
observed a significant deviation in the FWHM of the BAT PSF. Tab. 1 shows the dimension 
in physical units (e.g. kpc) of the extraction radius of 10' at the redshift of the source. 

The level of the background was evaluated using blank-sky observations^ (e.g. see Lumb 
et al. 2002; Read & Ponman 2003; Nevalainen et al. 2005) which are described in details 
in Carter & Read (2007). The black fields were selected from the same sky region as the 
observation under analysis and with a similar foreground absorbing column density. These 
fields were then reprojected to the source (sky) coordinates system and processed in a similar 



'^Blank-observations are described and made available at: http:/ /xmm. vilspa.esa.es/external/xmm_sw_cal/background/bla 
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way as the observation under analysis. For each cluster, a background spectrum has been 
extracted from the exact same 10' region as the cluster. To allow for different intensities 
of the background components (between the source and the 'background' observations) we 
renormalized the background spectrum by the ratio of the total emission in the 10-12 keV 
band in the annulus with inner and outer radii of 10' and 12' respectively (see also Molendi 
& Gastaldello 2009). 

The results of the background subtraction change slightly by varying this renormaliza- 
tion constant within its error, as well as changing the extraction annulus of the background 
or the thresholds for the removal of the flaring episodes (in the hard and the soft bands). We 
noticed that the spectral results are robust for variation of the aforementioned parameters 
if a systematic uncertainty of ~2 % is applied to our background subtracted spectra. We 
thus employed this systematic uncertainty connected to the background subtraction, in the 
0.5-8 keV band, when fitting the XMM- Newton data. Finally, all the spectra were rebinned 
in order to have a minimum of 50 counts {>7 a) per bin. 

The most distant clusters in our sample provide a test bed for checking the goodness 
of the background subtraction employed in this work. Indeed, for the Bullet cluster, PKS 
0745-19 and Abell 2390 is possible to find spatial regions of the EPIC-PN CCD that are the 

least contaminated by the cluster emission. These regions provide a clean way to determine 
the background spectrum which suffers from different systematic uncertainties^ with respect 
to the use of blank flelds. In all these cases we found that within the aforementioned 
systematic uncertainty, the two background subtracted spectra (e.g. the one that uses blank 
field observations and the one that uses part of the CCD which is not contaminated by the 
cluster emission) are in good agreement. One such example is reported in Fig. 2, which 
shows the two background subtracted spectra (generated with the two different techniques 
described above) of PKS 0745-19 fitted with a single temperature thermal model. As it is 
apparent from this figure, there is very good agreement between the two spectra. 



2.4. Point-like Sources in the XMM-Newton Fields 

Given the extent of its PSF, BAT is unable to discriminate the cluster's emission from 
that of nearby point sources which fall in the cluster's region. Thus, if present these sources 
would contaminate with their signals the cluster's emission as seen by BAT. For this reason, 
when analyzing XMM-Newton data jointly with BAT, we do not filter out the point sources 



^The main systematic uncertainty is that the background spectrum is extracted from a region different 
from the region used to extract the cluster's spectrum. 
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Fig. 2. — Spectra of PKS 0745-19 using two different background subtraction techniques: 
blank field observations (in red) and an extraction region free of cluster's emission (in black). 
The spectra were rebinned differently to simplify the comparison. 

which are clearly resolved and detected by XMM- Newton. However, we employ two different 
approaches to determine whether these sources contaminate in any way the total X-ray signal 
in the BAT band. 

In the more general approach we evaluate the contamination from point sources using the 
2XMM catalog (Watson et al. 2009) which contains all point sources detected in the XMM- 
Newton fields including the observations used in this work (e.g. 1). From the 2XMM catalog 
we select all the point-sources detected in the 0.2-12 keV band, with a likelihood probability 
of being spurious -ln(P)>15 and within 10' from the Swift-BAT centroid. The count rate 
of all the sources in the field of view has been summed and converted by extrapolation into 
15-55 keV fluxes in units of erg cm~^ s— 1. In this conversion we assumed that the spectra 
of these sources were represented by a power-law with a photon index of r=2, absorbed 
by a cold material with column- density Nh equal to the Galactic value. We believe that 
this approach will yield conservative estimates for the contamination to the 15-55 keV signal 
produced by point-like sources. Indeed, by assuming a power-law spectrum we implicitly 
assumed that all sources detected by XMM- Newton are AGN, while it is known that a large 
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fraction of the sources detected in clusters' region show a thermal spectrum (see Fig. 10 
in Finoguenov et al. 2004, and references therein) and thus will have a negligible >10keV 
emission. As shown in Burlon et al. (2010), the average broad-band 1-200 keV intrinsic^ 
spectra of AGN is compatible with a power-law with a photon-index of ~2.0. Our estimates 
are reported in Tab. 2. For most of the clusters in this analysis, the estimated contamination 
from point-sources is a factor ~ 50 below the total (cluster plus sources) emission measured 
by BAT (see Tab. 3 for details), and thus negligible. For only three clusters, PKS 0745-19, 
Abell 1795, and Abell 2390 the contamination might be relevant (still a factor ~ 5 below 
the total flux) . 

In the second approach, we study the XMM-Newton observations closely and we extract 
the spectrum of the brightest sources (up to 5) in the field. The characterization of their 
spectra in the 0.5-8.0 keV band allows us to make a solid prediction (i.e. without assump- 
tions) of the contaminating signal in the BAT band. The findings will be discussed case by 
case in the next sections. 



2.5. Joint Analysis of XMM-Newton and Swift/B AT Data 

Spectral analysis of XMM-Newton and Swift/BAT data has been performed, for all the 
clusters, using XSPEC (version 12.5.1 in Arnaud 1996). Observations of the Crab Nebula 
showed that, in principle, the inter-calibration of the two instruments is good within ~ 5%. 
Indeed, as reported in Kirsch et al. (2005), the Crab Nebula 0.3-10 keV spectrum as observed 
with the EPIC-PN can be modeled as an absorbed power-law with a photon index of 2.125, 
a normalization of 8.86 and an absorbing column density of 4.08 x 10^^ atoms cm~^. For 
BAT'', the Crab Nebula 15-200 keV spectrum is compatible with a power-law with a photon 
index of 2.15 and normalization of 10.17. Thus, for the 15-55 keV band (where most of 
the clusters' signal is concentrated for BAT) the two Crab Nebula parametrizations yield 
a flux of 1.21x10-^ erg cm^^ g-^ and 1.28x10"^ erg cm^^ g-^ for XMM-Newton and BAT 
respectively. Thus, we expect the inter-calibration of the two instruments to be close to 
unity (within ~5%). However, one must take into account that the Crab Nebula is a very 
bright target for XMM-Newton and in order to avoid pile-up problems the observations 
presented in Kirsch et al. (2005) were performed in "Burst Mode" (see Kirsch et al. 2005, 



^Photoelectric absorption caused by the circumnuclear material around the source will make the AGN 
spectra look harder in the < lOkeV band, but the > lOkeV continuum will be mostly unaffected (see Burlon 
et al. 2010, for details). 

^See http://swift.gsfc.nasa.gov/docs/swift/analysis/bat_digest.html for details. 
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Table 1. XMM- Newton observations of BAT clusters of galaxies 



NAME 


R.A.*" 


Decl.'" 


Date 


OBSID 


iLxposure 


Radius'^ 




(J2000) 


(J2000) 






(ks) 


(kpc) 


AbcU 85 


10.4303 


-9.3483 


2002-01-07 


0065140101 


10.0 


635.4 


Abell 401 


44.7395 


13.5837 


2002-02-04 


0112260301 


13.6 


823.8 


Bullet 


104.6176 


-55.8974 


2000-10-21 


0112980201 


46.7 


2460.0 


PKS 0745-19 


116.8758 


-19.3462 


2000-10-31 


0105870101 


28.3 


1100.0 


Abell 1795 


207.1856 


26.5928 


2000-06-26 


0097820101 


66.5 


702.6 


Abell 1914 


216.5013 


37.8071 


2002-12-18 


0112230201 


25.8 


1663.8 


Abell 2256 


256.0720 


78.6301 


2006-08-04 


0401610101 


50.4 


662.4 


Abell 3627 


243.6066 


-60.8348 


2004-09-19 


0204250101 


22.6 


204.1 


Abell 3667 


302.9667 


-56.8407 


2004-05-03 


0206850101 


67.3 


636.1 


Abell 2390 


328.4471 


17.7516 


2001-06-19 


0111270101 


23.1 


2079.6 



'^Swift/BAT coordinates are from the work of Cusumano et al. (2010). 

''Nominal XMM-Ncwton exposure before data screening. 

■^Extraction radius in physical units corresponding to the region of radius 10' around 
the BAT position used to extract the XMM-Newton spectrum of the cluster. 



Table 2. Summed X-ray emission from all the point-sources detected by XMM-Newton in 

the fields of the clusters. 



NAME 


Count-rate** 


# Sources 


Fluxo.2-12keV 


Fluxi5_55 keV 




count s~^ 




erg cm~^ s~^ 


erg cm-2 3-1 


AbcU 85 


0.23 


10 


5.61x10^13 


1.78x10-13 


AbeU 401 


0.10 


4 


3.72x10-13 


1.18x10-13 


Bullet 


0.10 


4 


2.38x10-1"^ 


7.55x10-" 


PKS 0745-19 


0.63 


14 


4.88x10-^2 


1.55x10-12 


Abell 1795 


1.76 


13 
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"^Summed count rate of all sources detected by XMM-Newton in the 0.2- 
12keV band. 
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for details) rather than the PuU-Prame Mode normally used for studying diffuse sources with 
XMM- Newton. Thus, the same inter-calibration between XMM-Newton and BAT might not 
necessarily apply in this case. However, Burlon et al. (2010) performed spectral fitting of 12 
faint AGN using XMM-Newton (with the EPIC-PN in Full-Frame Mode) and Swzft/BAT. 
In all these cases the intcr-calibration between the two instruments has been found to be 
compatible with unity. As a strategy in the spectral fitting presented in the next section, 
we employed a normalization constant to take into account differences in the calibrations of 
the XMM-Newton and BAT instruments. This constant has been fixed to 0.95 to take into 
account the different Crab Nebula spectra as observed with the two instruments. However 
we performed some tests with the joint datasets of this paper and found out that changing 
the intcr-calibration constant by ±10 % produces a negligible change (e.g. less than 1 %) in 
the best-fit temperatures and their uncertainties. Thus the results that will be presented are 
robust against variation of the aforementioned inter-calibration constant. The reason for the 
small variation of the temperature with the inter-calibration constant lies in the different 
signal-to-noise ratios of the two datasets. Indeed, in a joint fit the best-fit temperature is 
entirely constrained by the signal in the XMM-Newton band and a small variation of the 
inter-calibration constant (e.g. moving the less significant BAT data up or down around the 
best fit) does not change the results. We also checked that leaving this constant free to vary 
did not produce any appreciable improvement in the fit (in terms of goodness of fit) for all 
the clusters presented in this work. 

We started fitting all the spectra with a single-temperature thermal model (APEC) 
with absorption fixed at the Galactic value. Only if the value of the x^/ dof was significantly 
greater than 1, we tried to add a second thermal model or a power law. In this case we chose 
the model which produced the best improvement in the fit (evaluated using the F-tcst) and 
the best residuals. In all spectral fits all the parameters are tied together within the two 
datasets (e.g. XMM-Newton and BAT). 

In order to test which is the maximum level of non-thermal emission which is allowed by 
our data, a power law has been added to the best-fit model of every cluster. The power-law 
index has been fixed to 2.0, which is a value generally accepted for the non-thermal hard 
X-ray component generated by IC of relativistic electrons off CMB photons (e.g. Reimer 
et al. 2004; Nevalainen et al. 2004). We then let the power-law normalization vary until 
the Ax^ increment was larger than 2.7(6.64). According to Avni (1976), this gives the 90% 
(99%) confidence level on the parameter of interest. This allows us to investigate the level 
of non-thermal fiux which is consistent with our data. In the next sections, the details of 
the spectral analysis of each single cluster are reported. 
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3. Spectral Analysis of Individual Clusters of Galetxies 

3.1. Abell 0085 

Abell 85, a galaxy cluster at z=0.0521, has been detected for the first time in the X-ray 
band by Ariel V, thus it is one of the first galaxy clusters ever detected in X-rays (Mitchell 
et al. 1979; McHardy 1978). A detailed analysis performed with Einstein (Jones & Forman 

1999) showed that the ICM temperature is in the 7-9 keV range. ROSAT observations 
(Prestwich et al. 1995) revealed that the temperature and surface brightness structures of 
A85 are not regular, implying that the cluster is dynamically disturbed. At the same time 
the high central gas density indicates the presence of a cool core. Ensslin et al. (1998) showed 
that A85 has recently experienced a major merging, as indicated by the presence of a radio 
rehc. The estimate of the magnetic field intensity is, for the radio relic, of the order of 
B~2.6//G (Ensshn et al. 1998). Slee et al. (2001), using high-resolution radio observation, 
determined that the flux and spectral index of the relic at 1.425 GHz are SR=40.9mJy and 
a —3. 

The combined XMM- Newton and Swift/BAT dataset when fit by a single thermal model 
produces a x^/dof — 738.8/622, leaving unsatisfactory residuals at high energy. We thus 
added an additional thermal model to the fit which produces a significant improvement in the 

fit (Ax^ ^ 140 for 3 additional parameters). The best fit {x'^/dof= 602.1/619) temperatures 
are 6.09^929 keV and 1.72^QQg keV, while the respective metallicities are: 0.33too3 and 
0.15^003. The low temperature component accounts for the cool core of the cluster. Our 
results are in good agreement with the one reported by Durret et al. (2005) using XMM- 
Newton data alone. 

We also tried to add a power law model to the single thermal model. The fit improves 
with respect to the single temperature thermal model with a Achi'^ = 59 for 2 additional 
degrees of freedom. The best-fit temperature becomes in this case 5.09^Q :[gkeV and the 
photon index of the power law 2.68^0^3. However, the Ax^ is noticeable larger when the 
sum of two thermal model is fit to the data and we consider this model to be the best 
representation of our dataset (parameters reported in Tab. 3). Both spectral fits are reported 
in Fig. 3. 

We derive a 99 % upper hmit on any non-thermal component of 2.51 x 10~^^ erg cm~^ 
in the 50-100 keV band. The upper limit on the non-thermal luminosity in the 20-80 keV 
band is 3.62x10^^ erg s~^. When converted to the cosmology used by Nevalainen et al. 
(2004) (Ho=50km s"^ Mpc-^ qo=0.5 and A = 0) this becomes 6.68x10^^ erg s"^ which is 
in agreement with the value of 10.7lg;3 ^ 10^^ erg s~^ reported in Nevalainen et al. (2004). 
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Fig. 3. — Left Panel: Spectrum of Abell 85 fitted with the sum of two thermal models. 
Right Panel: Spectrum of Abell 85 fitted with the sum of a thermal model and a power 
law. 

3.2. Abell 401 

Abell 401, at z=0.074, is part of a cluster pair with Abell 399 which is in a pre-merging 
state (e.g. Karachentsev & Kopylov 1980; Fujita et al. 1996, and references therein). It is a 
rich cluster with a temperature of the ICM in the 7-8 keV range (Fujita ct al. 1996). Recently, 
using XMM-Newton, Sakelliou & Ponman (2004) found an average ICM temperature of 
7.23+ 

a2i keV. Abell 401 was one of the first clusters, along with Coma, that were discovered 
to host an extended radio emission (Harris & Romanishin 1974). This radio halo was recently 
confirmed by deep VLA observations (Bacchi et al. 2003). The intensity at 1.4 GHz was found 
to be = 17 ± 1 mJy 

The BAT spectrum is well fit by a brcmsstrahlung model with a plasma temper- 
ature of 7.79^2:86 keV. The combined XMM-Newton-BAT dataset is reasonably well fit 
(x^/(io/=766. 3/655) by a single thermal model with a temperature of 7.19±0.17keV and 
an abundance of 0.25±0.03 solar. Still, adding a second thermal model improves the fit 
substantially (Ax^ = 33.9 for 3 additional parameters). The temperature and abundance of 
the hot component are respectively 8.6llo;46 

keV and 0.301^^. Those of the cold component 
are respectively 2.05lo;45keV and 0.161q;J8. This fit is shown in the left panel of Fig. 4. 

We also tried a fit with the sum of a thermal model and a power-law (see right panel of 
Fig. 4). The best- fit temperature and power-law photon index are respectively 7.55lQ;27keV 
and 2.161^:^^. The improvement in the goodness of fit, with respect to the single thermal 
model, is Ax^ =9.1 for 2 additional parameters and is clearly not a better fit than the sum 
of two thermal models. Indeed the F-test yields a probability of 0.02 and 1.7x10"^ for the 
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power-law and the additional thermal model respectively of being spurious. For this reason 
wc believe that the sum of two thermal model is a more adequate representation of the 
XMM-Newton/BAT dataset and we report its best-fit parameters in Tab. 3. 

Since no spectral index for the radio emission is available in the literature we adopt a 
value of a =2.0. Using a power-law with a photon index of 2.0 we derive that the 99% CL 
upper limit on the non-thermal component in the 50-100 keV band is 2.2x10"^^ erg cm~^ 
s~^. As a final note, the brightest point-source in the XMM-Newton 10' region is located at 
R.A.(J2000) = 02:59:05.5 Decl.(J2000) = 13:39:44.9. Its spectrum is very soft and consistent 
with a bremsstrahlung model with a temperature of 0.3±0.1keV. Its fiux in the 2-10 keV 
band is 1.9x10"^® erg cm~^ s^^ and thus it is negligible when compared to the cluster signal 
in both the XMM-Newton and BAT bands. 




Fig. 4. — Left Panel: Spectrum of Abell 401 fitted with the sum of two thermal models. 
Right Panel: Spectrum of Abell 401 fitted with the sum of a thermal model and a power 
law. 



3.3. Bullet cluster 

IE 0657-56 is a distant cluster {z = 0.296), originally detected in the Einstein survey 
(Tucker et al. 1995). ROSAT and ASCA have shown that the Bullet is one of the hottest 
(kT=17.4±2.5 keV), and most massive cluster known (Tucker et al. 1998). The same data 
show that IE 0657-56 is undergoing a major merger process. Liang et al. (2000a) found out 
that IE 0657-56 contains a very luminous radio halo whose surface brightness closely follows 
the X-ray one. Weak and strong lensing reconstruction of the Bullet Cluster are one of the 
best evidence for the existence of dark matter (e.g. Clowe et al. 2004; Markevitch et al. 2004; 
Clowe et al. 2006; Bradac et al. 2006). The Chandmhigh resolution image of a bullet-like gas 
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cloud moving in the cluster core with a bow shock front, gained IE 0657-56 the name Bullet 
Cluster (Markevitch et al. 2002). The average temperature they report ranges from 14-15 
keV to more than 20 keV. Deeper Chandra observations showed that, away from the bullet, 
the radio halo peak is offset from the X-ray peak, which is centered on region hosting the 
hottest gas (Govoni et al. 2004). The Bullet Cluster was observed also with XMM-Newton 
(Zhang et al. 2004; Finoguenov et al. 2005; Zhang et al. 2006) and RXTE (Petrosian et al. 
2006). The latter determined that the spectrum of the Bullet cluster can be fit equally well 
by the sum of two thermal models or by the sum of a thermal and a power-law model. They 
also estimated that the equipartition value of the magnetic field intensity is ~1.2//G. IE 
0657-56 has a complex radio morphology. The diffuse radio halo detected by Liang et al. 
(2000a) has a flux density Sr =78 mJ at 1.3 GHz with a spectral index is a = 1.2. 

The combined XMM-Newton and BAT data can be successfully fit (xV^^o/=524.7/513) 
be a single-temperature thermal (APEC) model (see left panel of Fig. 5). The best-fit 
temperature is 12.57lQ;g5keV while the abundance is 0.25^0 Qg solar. The temperature is 
in moderate good agreement with the values of 14.5lfykeV and 14.8lj^'2keV as observed 
by ASCA and Chandra respectively (Liang et al. 2000b; Markevitch et al. 2002) and is 
well contained in the range of temperatures observed with Chandra (see above references). 
Following Petrosian et al. (2006) we added a power-law model to the fit. The fit improves 
(Ax^ Ri23 for 2 additional parameters, i.e. 4.4 cr) and the best-fit temperature becomes 
14.77^0^72 keV (in agreement with ASCA and Chandra results), while the power- law index 
is 1.861q;i4 ■ This model fit is reported in the right panel of Fig. 5. The non-thermal 
20-100 keV fiux is ?>At\l x 10-^^ 

erg cm s . These values are in good agreement with 
those reported by Petrosian et al. (2006). Moreover, recently Million & Allen (2009), using 
Chandra, reported the detection of non-thermal ffux in the 0.6-7.0 keV band at a level of 
0.95+ 

an ^ 10 erg cm ^ s ^. From our analysis we derive that the 0.6-7.0 kcV power-law 
fiux is 2,.2,t\l X 10-^2 

erg cm ^ s ^ and thus slightly brighter than their reported fiux. 

For the sake of completeness we also tried to fit the spectrum of the Bullet cluster 
with the sum of two thermal models. The hot, most intense component, shows as before 
a temperature of 15.41^5 keV and an abundance of O.SOlo'gg solar. The second component 
displays a temperature of l.llQ;2keV, an abundance compatible with zero and a 1-10 keV 
fiux of 6.3j^2 i ^ 10~^^erg cm~^ s~^. This fit is shown in the lower panel of Fig. 5. Both 
fits represents a reasonable description of the data. The thermal plus power-law model is 
slightly worst (x7^^o/=501.7/511) than the sum of two thermal models (xV^^o/=499.9/510). 
However inspection of Fig. 5 shows that the thermal plus power-law model explains better 
the residuals at high-energy, albeit the BAT data are not very significant above 50 keV. We 
also checked that leaving the intcr-calibration between BAT and XMM-Newton free to vary 
(see § 2.5) does not changes the results presented here. 
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Finally we also verified the contribution of point sources to the overall signal. We found 
that only two sources produce a signal comparable to the excesses seen here. The first one is a 
bright point sources locate south-west of the cluster core at a position R.A.(J2000)=06:58:13.8 
and Decl.(J2000)=-55:59:20.6. Its spectrum can be fit by an absorbed bremsstrahlung model, 
where the absorption is compatible with the Galactic one and the temperature of the plasma 
is S.Ql^l'alkeV. The 1-10 keV flux is 2.90t\il x 10"^^ erg cm'^ s'^ The second one is lo- 
cated at R.A.(J2000)=06:58:03.8 and Decl.(J2000)=-56:01:13.1 and its spectrum can be fit 
with an absorbed power law where the absorption is in excess of the Galactic one with 
Nh =9.3l4:I X lO^^cm-^ and the photon index is 1.46lo:3i- When extrapolated to the 20- 
100 keV band the source flux is 5.0!t2!45 ^ 10~^^erg cm~^ s~^. It is thus clear that both 
sources cannot account for the observed signals, indeed their fluxes are a factor ~10 below 
the flux of the 'cold' component and the non-thermal component seen in the spectral flts 
described above. 

Moreover, we checked whether the results reported above might be connected to some 
residual background contamination in XMM-Newton which was not accounted for correctly 
in the extraction of the background spectrum from blank fleld observations. For this purpose, 
we extracted a background spectrum from a region of the XMM-Newton CCD which is the 
least contaminated by the cluster emission. Since the Bullet cluster is at moderately high- 
redshift this is possible. We also extracted the background spectrum from a region whose 
area was the same as that one used to extract the cluster emission and whose position is 
diametrally opposed to the cluster region with respect to the pointing direction (to ensure a 
similar effective area over the two regions) . The analysis of the Bullet cluster spectrum using 
this background strategy conflrms the above results. In particular both the 'cold' component 
and the non-thermal power law are confirmed and the derived fiuxes are consistent with those 
reported above. Thus, we can exclude an instrumental origin for both components. 

In order to understand whether both components co-exist, we tried a fit with a model 
which is the sum of two thermal component and a power law. Not surprisingly this model 
produces a good fit to the data (x^/nc(/'=498. 1/508). All parameters, with the exception of 
the normalization of the 'cold' component, are fully compatible and in good agreement with 
the parameters reported above. Indeed, from this best-fit we derive that the cold component 
has now a 1-10 keV flux of 7.7^.]'} ± 10~^^erg cm~^ s~^ and thus a factor of ~ 10 below 
the flux derived from the flt using two thermal models only. Moreover, as the statistical 
uncertainty shows, this component is now compatible with zero at 90 % confldence. 

On the other hand, the parameters of the power-law component are robust with respect 
to the variation of the other parameters (see Fig. 6 ). For this and all the other reasons 
explained above we believe that the description of the cluster spectrum in terms of a single 
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thermal model and a power law is the best and most reliable one. The parameters of this fit 
are reported in Tab. 3. Our analysis thus confirms the presence of a power-law component 
in the spectrum of the Bullet cluster as reported by Petrosian et al. (2006). 
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Fig. 5. — XMM-Newton and BAT data for the Bullet cluster fitted with: 1) a single thermal 
model (upper left), 2) the sum of a thermal and a power-law model (upper right), 3) the 
sum of two thermal models (bottom left), and 4) sum of two thermal models (thin and 
long-dashed lines) and a power law (bottom right). 



3.4. PKS 0745-19 

Early Einstein and EXOSAT observations showed that PKS 0745-19 (also known as 
4U 0739-19) is one of the largest cool core cluster known (Fabian et al. 1985; Arnaud et al. 
1987; Edge et al. 1990; White et al. 1997). The ICM temperature as measured with ASCA 
and ROSAT (Allen et al. 1996; Pierre & Starck 1998; Peres et al. 1998) agrees well with 
that one found with BeppoSAX, of about 8.3 keV (De Grandi & Molendi 1999, 2001, 2002). 
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Fig. 6. — Contour (1,2, and 3 a) plots of the parameters of the power-law component versus 
the parameters of the most intense thermal component for the Bullet cluster. 

Chen et al. (2003), using XMM- Newton and Chandra, reported an enhanced diffuse X-ray 
emission in correspondence of bright radio lobes. They discuss the possibility of buoyant 
bubbles to explain the observed X-ray and radio emission. However Dunn & Fabian (2006), 
using Chandra data, find "no clear" evidence of such radio bubbles. Ball et al. (1993) and 
Baum & O'Dea (1991) classify PKS 0745-19 as an amorphous radio source, displaying both 
a compact radio source (five times brighter than Perseus) and a diffuse emission. These 
unusual radio properties do not seem just the result of AGN activity, but the result of a 
merger and/or buoyant plumes. Baum & O'Dea (1991) report that the intensity of the 
diffuse flux at 5 GHz is 265 mJy, while the spectral index is a = —1.4. For these values, the 
equipartition magnetic field is in the 20-50 fiG range. 

Our data point to a thermal origin of the hard X-ray emission. Indeed, the XMM- 
Newton and BAT spectra can be successfully modeled by a single temperature thermal model 
(x^/c?o/=610. 5/581) with a temperature of 6.69]'^Q;27keV and an abundance of 0.35^;g;o3. 
However, adding a second thermal component produces a noticeable improvement (Ax^ = 
23.2 for 3 additional degrees of freedom, corresponding to 5.0x10"^ chance). The hot most 
intense component has now a temperature of 7.96^o;54keV and an abundance of 0.40lo;o5, 
while the cold component displays a temperature of 2.16;^^o!56 ^^V and an abundance of 
0.3llo;oQ. This fit is shown in the left panel of Fig. 7. The temperatures observed here are 
consistent with those found by George et al. (2009) using Suzaku. 

We also tried a fit with the sum of a thermal and a power-law model. The temperature 
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and abundance of the thermal component are 7.531° ggkeV and an abundance of O-STIqIo?) 
while the photon index of the power law is 4.39]'^2 qi- The improvement in the with 
respect to the single temperature thermal model is Ax^ = 17 for 2 additional parameters 
(corresponding to 0.2 % chance). This fit is reported in the right panel of Fig. 7. Nevertheless, 
given the very soft spectral index (which accounts for the 'cold' component) and the marginal 
improvement in the fit statistics, we believe that the sum of two thermal models is a more 
rehable interpretation of the XMM-Newton/BAT dataset. The parameters of this fit are 
thus summarized in Tab. 3. 

Using the sum of two thermal models as a baseline spectral fit we estimated the 99 % 
CL to a non-thermal component in the 50-100 keV band using a power law with a photon 
index of 2.0. This is found to be 1.55 xlO~^^ erg cm"^ s""^. If instead a photon index of 2.4 is 
used, the upper limit in the 50-100 keV band becomes 1.10x10"^^ erg cm~^ s""^. Finally, the 
brightest point-source in the XMM-Newton 10' region is located at R.A.(J2000) = 07:47:19.0 
Decl.(J2000) = -19:24:02.3. Its spectrum is very soft and consistent with a bremsstrahlung 
model with a temperature of 0.7±0.3keV. Its fiux in the 2-10 keV band is 2.6x10"^^ erg cm~^ 
s~^ and thus it is negligible when compared to the cluster signal in both the XMM-Newton 
and BAT bands. 




Fig. 7. — Left Panel: Spectrum of PKS 0745-19 fitted with the sum of two thermal models. 
Right Panel: Spectrum of PKS 0745-19 fitted with the sum of a thermal model and a 
power law (dashed line). 
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3.5. Abell 1795 

Abell 1795 is a compact and rich cluster with a strong cool core (e.g., Edge et al. 1992; 
Briel & Henry 1996; Markevitch et al. 1998; Tamura et al. 2001). It has been extensively 
observed with HEAO-1, Emstem, EXOSAT, ROSAT, BeppoSAX, XMM-Newton, Chandra 
and RXTE (e.g. Kowalski et al. 1984; Rhee & Latour 1991; Edge et al. 1990; Arnaud et al. 
1991; Markevitch et al. 1998; Sanders et al. 2000; Arnaud et al. 2001; Nevalainen et al. 2004; 
Revnivtsev et al. 2004; Vikhlinin et al. 2005, and references therein). EGRET was also used 
to set an upper limit to the 7-ray emission from the ICM (above 100 MeV, Reimer et al. 
2003). Outside the cool region, the average temperature is of 6-7 keV (e.g. Arnaud et al. 
2001). Abell 1795 is a relaxed cluster, although there is evidence that the central brightest 
galaxy is not at rest (Hill et al. 1988) and that there is inner gas sloshing in the potential well 
(Markevitch et al. 2001; Ettori et al. 2002). No hard X-ray excess has ever been reported 
for this cluster. Fabian et al. (2001) discovered a 40" long cold filament in the core of the 
cluster, in coincidence with an Ha filament. The straightness of the filament indicates that 
the ICM is not very turbulent. Abell 1795 has two X-ray dim regions (e.g., Fabian et al. 
2001; Ettori et al. 2002) in correspondence of radio bubbles (Dunn et al. 2005). The radio 
morphology is dominated by two radio regions (Guthrie 1974; Owen 1975; Dagkesamanskii 
et al. 1982; Alikberov et al. 1983; Burns 1990; Owen et al. 1993; Owen & Ledlow 1997; Dunn 
et al. 2005), but there is no strong evidence for a large scale radio halo (Hanisch 1982). 
Ge & Owen (1993) report Faraday rotation measurements of the small central radio galaxy 
1346+268. They conclude that the associated magnetic field must be >20 /xG and that it is 
most likely associated with the ICM rather than with the small radio source. 

Using the parameters of the surface brightness reported by Briel & Henry (1996) (e.g. 
core radius of 5.15' and (3 =0.93), we derive that our standard selection region of 10' radius (in 
XMM-Newton) includes 97-99% of the cluster emission. The XMM-Newton and Swift/BAT 
data are best modeled {x^/dof — 892.1/1275) by a single thermal model with a temperature 
of 4.82lQ;j^5 keV and an abundance of 0.451q;o4. Our results are in good agreement with 
the ones of Tamura et al. (2001). Given the very good x^, adding other models (or free 
parameters) does not improve the fit. Thus, we believe that the single thermal model is a 
good representation of the XMM-Newton/BAT dataset. This fit is shown in the left panel 
of Fig. 8. 

As a final note, the field of Abell 1795 contains a bright AGN located at R.A.(J2000)=13:48:35.2 

and Decl.(J2000)=26:31:08.9, which lies less than 3' from the cluster core. This source is 
the Seyfert 1 galaxy IE 1346+26.7. In XMM-Newton it exhibits an unabsorbed power-law 
spectrum with a photon index of 2.34+0.04 and a 2-10 keV flux of 8.7(+0.1) x 10^^^ erg 
cm~^ s~^ (approximately 1 % of the cluster's signal in that band). When extrapolated to 
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the 15-55 keV band, the AGN flux becomes 3.7(±0.1) x 10~^^erg cm-^ and thus still 
a factor ~ 5 below the cluster emission in that band (see Tab. 3). We checked that the 
results reported above do not change if the AGN emission is properly modeled (e.g. with the 
parameters of the AGN power-law allowed to vary within their errors) in the joint fit to the 
overall cluster's spectrum. However, when computing the (50-100 kcV) 99% CL upper limit 
on the non-thermal component (using a power-law with a photon index of 2.0) we derive 
that this is 1.57x10"^^ erg cm~^ s~^ if the AGN contribution is not taken into account, or 
1.06x10"^^ erg cm~^ s~^ if it is. We will thus use this second upper limit in Tab. 4. 




Fig. 8. — Left Panel: Spectrum of Abell 1795 fitted with a single temperature thermal 
model. Right Panel: Spectrum of Abell 1914 fitted with a single temperature thermal 
model. 



3.6. Abell 1914 

Abell 1914 is a regular and smooth galaxy cluster. It has been observed with ROSAT 
(Ebcling ct al. 1996; Buote & Tsai 1996; Bohringer et al. 2000), ASCA (White 2000; Ikebe 
et al. 2002) and Chandra (Govoni et al. 2004; Baldi et al. 2007). By means of a comparison 
of X-ray and radio maps, Govoni et al. (2004) discuss a possible merger scenario. Baldi et al. 
(2007) report an average ICM temperature of 9.20 ± 0.39 keV. Abell 1914 is known to host 
a very steep radio source (Kulkarni et al. 1990) and a radio halo (Giovannini et al. 1999; 
Kempner & Sarazin 2001; Bacchi et al. 2003). The point sources make the estimate of the 
diffuse flux density difficult. We adopt the value of S'j?=64mJ at 1.4 GHz and the reported 
spectral index of a = 1.8 from Bacchi et al. (2003). The equipartition magnetic field is 
0.5 A^G. 

The Swift/ BAT spectrum can be fit by a bremsstrahlung model with a plasma tempera- 
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ture of IMttll keV. The combined XMM-Newton - Swift/BAT dataset can be successfully 
modeled (x^/'^o/=355. 1/351) with a single-temperature thermal model (see right panel of 
Fig. 8). The best-fit temperature and metallicity are 11.14l}'Jg keV and 0.19±0.14 solar, in 
agreement with the studies mentioned above. Given the good adding other models to 
the single thermal model does not improve the fit results. The 99 % CL upper limit on the 
50-100 keV non-thermal flux, evaluated with a power law with a photon index of 2.0, flux 
is 1.08x10"^^ erg cm~^ s~^. If we use a power law with a photon index of 2.8 (in line with 
the radio photon index) then the upper limit is much tighter and it becomes 4.60x10"^^ erg 
cm~^ s~^. However, this upper hmit to the non-thermal flux in the (BAT) 50-100 keV band, 
is entirely driven by the XMM-Newton signal below 2 keV. Indeed, if we repeat the same 
process described above, but using only Swift/BAT data, then the 99 % CL upper limit (us- 
ing a photon index of 2.8) is 1.16x10^^^ erg cm~^ s^^ and thus in line with the one computed 
using a power-law with a photon index of 2.0 and the entire XMM-Newton/BAT dataset. 
We thus believe that this (e.g. 1.08x10"^^ erg cm~^ s~^) is for the 50-100 keV band a more 
reliable upper limit. 



3.7. Abell 2256 

Abell 2256 is a rich cluster at a redshift of 0.0581, bright both at radio and X-ray energies 
(e.g. Bridle & Fomalont 1976; Briel et al. 1991; Henriksen 1999). It has been studied several 
times at X-rays and the disturbed morphology of the X-ray temperature map indicates a 
cluster in an advanced merging stage (e.g. Molendi et al. 2000). Abell 2256 is one of those 
clusters for which a claim of significant detection of non-thermal emission has been reported. 
Indeed, Fusco-Femiano et al. (2000), using data from BeppoSAX, reported the detection of a 
hard X-ray excess at the 4.6 cr level. The 20-80 keV flux of this excess is 1.2x10"^^ erg cm~^ 
s"-*^. Rephaeli & Gruber (2003), using RXTE data, reported the detection of an hard X-ray 
excess whose 20-80 kcV flux is 4.3^*^40 x 10~^^erg cm~^ s~^ (errors are 90% CL) and thus a 
factor of ~3 fainter than the one reported by Fusco-Femiano et al. (2000), but marginally 
consistent with it. A re-analysis by Fusco-Femiano et al. (2005) conflrmed the BeppoSAX 
detection (at 4.8 a) albeit at a lower 20-80 keV flux of 8.9+3:g x lO'^^erg cm'^ s'^ In the 
radio band, Abell 2256 displays an extremely complex morphology consisting of a bright 
rehc and a fainter steep-spectrum radio halo located in the cluster center (Clarke & Ensslin 
2006). The total flux density of this radio halo is 100 mJy at 610 MHz (Rengelink et al. 1997) 
and 103.4 mJy at 1369 MHz (Clarke & Ensslin 2006), while the spectral index is q;=1.8. 

The BAT data alone are well flt by a bremsstrahlung model with a temperature of 
9.8l3;8 keV. Adopting the values reported by Briel et al. (1991) (e.g. core radius of 4.83±0.17 arcmin 
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and ^=0.756±0.013), we derive that selecting photons within 10' of the core includes ~95% 
of the cluster's emission. The joint XMM-Newton-BAT datasct is well fit {x^ / do f =445. 6/ 445) 
by a single thermal model with a temperature of 8.84tQ g5 keV and an abundance of 0.22±0.06. 
The best fit is shown in Fig. 9 while the parameters are reported in Tab. 3. Henriksen (1999) 
found out that the best spectral model reproducing the RXTE/ASCA datasets, for Abell 
2256, is produced by the sum of two thermal models. In that work the hot and the cold 
components have a temperature of ~7keV and ~lkeV respectively. Following his exam- 
ple, we added a second thermal model to the fit keeping the abundance of this additional 
component fixed at 0.3 (allowing this parameter to vary does not change the results). The 
best-fit temperature of the additional component is 1.08±0.39keV in good agreement with 
the results of Henriksen (1999) while the temperature and abundance of the hot component 
did not vary appreciably. However, the improvement in the Ax^ is 2.5 for two additional 
parameters and thus not significant (i.e. the probability that the improvement was obtained 
by chance is ~0.3). Wc thus believe that the single temperature thermal model discussed 
above represents the best description of the XMM-Newton-BAT dataset. 

The 99 % CL upper limit on the 20-80 keV non-thermal fiux is 6.1x10"^^ erg cm~^ s~^. 
The upper limit derived by our analysis is lower than the hard X-ray excess claimed by 
Fusco-Femiano et al. (2005). Even using BAT data alone, the 99% upper limit in the 20- 
80keV band is 4.6x10"^^ erg cm~^ s~^ and thus inconsistent with the BcppoSAX result (but 
not with the RXTE one). In our band (50-100 keV) the upper hmit, derived from the joint 
dataset, is 2.41 xlO~^^ erg cm~^ s~^. If instead of a power law with an index of 2.0, we use 
a power law with a photon index of 2.8 (the value of the radio halo) the 50-100 keV upper 
limit on the non-thermal emission would be 1.97x10"^^ erg cm~^ s~^. 



3.8. Abell 3627 

Abell 3627, known also as Norma cluster at z=0.015, is a nearby massive cluster located 
behind the Milky Way in the core of the Great Attractor and discovered as an important 
component of the local large-scale structure by Kraan-Korteweg et al. (1996). It is a very rich 
cluster with a mass comparable to that of Coma and Perseus (i.e. >2xl0^^ M©). Early X- 
ray observations with ROSAT and ASCA show that the cluster is not spherically symmetric 
and has a strong temperature gradient (A/cT ~3keV) in the direction of the elongation 
(Boehringcr et al. 1996; Tamura et al. 1998). This fact indicates that Abell 3627 is in the 
stage of a major merger. This cluster also exhibits spectacular head-tail radio-galaxies (e.g. 
Sun et al. 2010) which are galaxies hkely traveling at high velocities through the ICM (see 
e.g. Sarazin 1988). At the cluster center, PKS 1610-608 displays, in radio, two powerful 
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Fig. 9. — Spectrum of Abell 2256 fitted with a single temperature thermal model. 

jets and two lobes whose surface brightness peaks respectively at ~1' and ~5' away from 
the galaxy. Jones & McAdam (1996) find that the intensity of the magnetic field, derived 
assuming equipartition, is ~15 /iG at the position of the jets and ~ 5 /iG at the lobes. 

Considering the values for the surface brightness reported by Boehringer et al. (1996) 
(e.g. core radius of 9.95' and /3 =0.55) and Fig. 1 it is clear that Norma should be detected, 
by BAT, as an extended source. However, its elongation (e.g. not being spherical symmetric) 
does not allows us to determine a-priori the expected hkely flux suppression in BAT. This 
is made even more complex by the presence of a nearby AGN (IGR J161 19-6036) which is 
detected at a significance of ~ lOc by BAT (Cusumano et al. 2010). This AGN is located 
at ~20' away from the BAT centroid of Norma and the two sources appear separated. 
Fig. 10 shows the contours of the surface brightness of Norma (as derived from ROSAT- 
PSPC observations) superimposed on the BAT significance map for that region. It is clear 
from the ROSAT contours that Norma extends hkely all the way to the nearest point-like 
source (IGR J16119-6036). We should thus expect a contamination of the cluster thermal 
emission in the BAT spectrum of IGR J161 19-6036. This would imply that BAT detects the 
Norma cluster as an extended source. We started fitting the BAT data alone for the Norma 
cluster with a bremsstrahlung model. The fit is acceptable (x^/dof—lAA/lA) and the best- 
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fit temperature is 11.6]'^3 3 keV. Next, we extracted Swift/XKT data for IGR J16119-6036 and 
fitted them together with the BAT data for this source. The resuhs are reported in Fig. 11. 
When using a single absorbed power law (which fits the XRT data alone well), the fit to the 
XRT-BAT dataset is unacceptable with x^/dof = 31.5/16 and leaves (as it can be seen in 
the left panel of Fig. 11) residuals in the BAT band. We then added a bremsstrahlung model 
(only for the BAT, since XRT detects IGR J16119-6036 as point-source) to the AGN power 
law, to check whether the BAT spectrum is contaminated by the thermal emission from the 
Norma cluster. The fit with this model is good {x^ /do f =12.1/14) and the photon index of 
the power law is 1.67±0.16 while the temperature of the thermal model is 12.3±4.2keV. The 
temperature is in good agreement with the temperature of the Norma cluster (as measured 
with BAT) reported above and thus we conclude that there is significant contamination, 
particularly below 50keV, of cluster's emission in the BAT spectrum of IGR J16119-6036. 
This also means that BAT detects the Norma cluster as an extended source. 
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Fig. 10. — Contours of the surface brightness of the Norma cluster, as derived from ROSAT- 
PSPC observations, superimposed on the Swift/BAT significance map. 
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Fig. 11.— Left Panel: XRT and BAT spectrum of IGR J16119-6036 fitted with a power 
law. Right Panel: XRT and BAT spectrum of IGR J16119-6036 fitted with the sum of a 
power law and a thermal model (the latter only for the BAT data). 



Given the finding that BAT 'resolves' Norma, it becomes difficult to determine how 
much fiux has been suppressed by the mask and such detailed analysis will be left to a 
future paper. For this reasons, instead of providing a joint fit to XMM-Newton and BAT, 
we perform two separate spectral fits. The BAT data, as already described, are shown in 
the left panel of Fig. 12 and are well fit with a bremsstrahlung model with a temperature 
of ll.ellgkeV. The XMM-Newton data (extracted around 10' from the BAT centroid) are 
well fit Ix^/ do f ^402.6/377) by an APEC model with a temperature of 5.53l^:|keV and 
an abundance of 0.26to'o3. This spectrum is shown in the right panel of Fig. 12. Clearly, 
this analysis points to a difference in the temperature of the plasma as measured with the 
two instruments. This piece of evidence^ would point towards the existence of regions of 
hot gas in the Norma cluster. While spatially-resolved spectroscopy is not available for this 
cluster, both Boehringer et al. (1996) and Tamura et al. (1998) find that for some regions of 
the cluster temperatures as high as 7-10 keV might exist, thus in agreement with the BAT 
detection. For this cluster, we report in Tab. 3 the parameters of the best fit to the BAT 
data alone. Since it is resolved by BAT, part of its > 15keV fiux is lost in the background 
and an upper limit to the non-thermal emission will not be computed. 



The XMM-Newton and BAT temperature are still compatible with each other within ~2 a. 
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Fig. 12.— Left Panel: BAT data for Abell 3627 fitted with a bremsstrahlung model. Right 
PanelrXMM-Newton data for Abell 3627 fitted with a single thermal model. 

3.9. Abell 3667 

Abell 3667 is a cluster at z=0.055 discovered by HEAO (Piccinotti et al. 1982). ROSAT 
revealed that Abell 3667 is a dynamically interacting system with a significant X-ray emis- 
sion associated with a group of galaxies which is hkely merging with the cluster. ROSAT 
measured for the cluster an average temperature of ~6.5keV (Knopp et al. 1996). Its dy- 
namically complex structure has been investigate by Vikhlinin et al. (2001) using Chandra, 
revealing that the dense cool core is moving with high velocity through the hotter, less dense, 
surrounding gas, creating a cold front. Vikhlinin et al. (2001) estimated the intensity of the 
magnetic field, in the vicinity of the shock region, to be B~10 /xG. The magnetic field near 
the cold front is expected to be stronger and to have a very different structure compared to 
the bulk of the ICM. These peculiar characteristics make A3667 a good candidate for the 
detection of a hard X-ray excess of non-thermal origin. This component has, indeed, been 
reported in hard X-ray spectrum measured by Beppo-SAX (Fusco-Femiano et al. 2001). 

Abell 3667 was recently studied in detail by Nakazawa et al. (2009) using data from 
Suzaku. When modehng the XIS (0.7-8.0 keV) and the HXD (15-40 keV) spectra they found 
that a single thermal model fails to explain the hard X-ray data and that another component 
is needed. This is required to be a very hot thermal component with T=19.2l4;QkeV or a 
power-law with an index of 1.391q;i7 (Nakazawa et al. 2009). 

Fig. 13 shows the BAT significance map with superimposed contours from ROSAT 
(X-rays) and SUMSS (radio). It is clear that the BAT detection is associated with the 
core of the cluster and it is not compatible as coming from the radio relic which hes >12' 
north-west of the BAT centroid. The BAT spectrum (reported in Fig. 14) shows that Abell 
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3667 is indeed an interesting cluster. A simple bremsstrahlung model fits the data well 
and the best fit temperature is IS.Tlg^g^kcV which is unusually high even for BAT and in 
agreement with the one found by Suzaku. Adopting the values for the surface brightness 
reported by Boehringer et al. (1996) (e.g. core radius of 2.97' and /3 =0.55) we derive that our 
standard selection, in XMM-Newton, of photons within 10' includes ^85-95 % of the cluster's 
emission. When analyzing jointly XMM-Newton Swift/BAT data we find that a single 
thermal model yields a best-fit temperature of 5.68±0.19 keV and an abundance of 0.21±0.04 
solar. This fit is reported in the upper panel of Fig. 15. It is apparent that this fit leaves 
unsatisfactory residuals at high-energy. We then tried adding a second thermal component. 
The fit improves and the two thermal components show a temperature of 13.5l2!2keV and 
3.9i2!ikeV respectively. The F-test shows that the probability of the second component to 
be spurious is only 4.86x10"^. 

An equally good fit can be obtained with the sum of a thermal and a power-law model. 
In this case the best-fit temperature is 5.91±0.05keV and the photon index is 1. 83^0 34- 
Again the F-test shows that the probability of the second component to be spurious is very 
low, 3.46x10"^ (e.g. the significance of the model is ~4.6(t). The power-law flux in the 10- 
40keV band is 4.91^° 00 ^ 10~^^ erg cm~^ s~^, a factor of ~10 fainter than the one reported 
by Nakazawa et al. (2009). 

The two models discussed here produce the same result in term of goodness of fit and 
point to the existence of a very hot region with a temperature ~ 13 kcV (as the BAT data 
alone testify). On a pure statistical basis the model with less parameters should be chosen 
(e.q. the thermal plus power law model). However, on a physical basis it is difficult to 
understand whether this excess is due to a hot component (as seen in other cases) or to a 
truly non-thermal power-law like one. We believe that the hot component is the more realistic 
hypothesis for several reasons: first the temperature of this component is not unusually high 
for massive and merging galaxy clusters. Second, Abell 3667 is known to have radio relics, 
but not a central radio halo (Rottgering et al. 1997). Since the BAT centroid (see Tab. 1) is 
compatible with the cluster core and not with the relic, this renders the interpretation of the 
hard X-ray excess as being non-thermal unlikely. In support of the thermal interpretation 
of the hard X-ray excess we note that the BAT centroid coincides with the hottest region 
(kT>8keV) as found with XMM-Newton (Briel et al. 2004). For this reason we decided to 
report the sum of the two thermal models in Tab. 3. 

In order to exclude that this hard X-ray excess originates from one of the point sources, 
we extracted the spectrum of the brightest X-ray sources located in the cluster field. Among 
all of them the brightest is the source positioned at RA(J2000)=303.14908 and Decl.(J2000)=- 
56.89704 with an uncertainty of 3" . The XMM-Newton spectrum is consistent with a sim- 
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pie power-law with an index of 1.73±0.20. Its flux extrapolated to the 10-40 keV band is 
cm ^ s ^ Since this flux is a factor > 10 fainter than the hard X-ray 
excess, we can exclude that the hard X-ray excess is caused point-like sources. Wc thus 
conclude, partly confirming the result of Nakazawa et al. (2009), that our data requires a 
hot component (kT=13.5^2 2 ^eV) or a power law with a photon index of I.83I034. In this 
last case the 50-100 keV non-thermal flux is 2.98^q;73 x 10~^^erg cm~^ s~^. 

Finally we note that, as can be seen from Fig. 13, there is a small ~3.5o' fluctuation 
in the BAT map ~12' North-East of the cluster core. However, we remark that until the 
5 a threshold is exceeded this has to be considered a statistical fluctuations. Indeed, the 
probability of observing a pure >3.5 a statistical fluctuations in the BAT map is quite large 
(i.e. 2.3x10"^) leading^ to a total of ~6700 statistical fluctuations. Moreover, no known 
AGN or (bright) X-ray sources are reported within 5' of this fluctuation and inspection of 
all the available X-ray data (XMM-Newton, Swift/XRT, ROSAT, etc.) did not reveal any 
potential candidate that might be the counterpart of this sub-threshold object. We thus 
believe this to be just a statistical fluctuation. 

3.10. Abell 2390 

Abell 2390 is a rich lensing galaxy cluster with a massive cool core (Pierre et al. 1996, 
e.g,). It is among the ten X-ray brightest galaxy clusters at redshift larger than 0.18 (Ebeling 
et al. 1996, e.g,). It has been observed with HEAO 1 and 2 (Johnson et al. 1983; Kowalski 
et al. 1984; Wood et al. 1984; Ulmer et al. 1986), Emstem (McMillan et al. 1989), ROSAT 
(Pierre et al. 1996; Ebeling et al. 1996; Pierre & Starck 1998; Rizza et al. 1998; Bohringer 
et al. 2000), ASCA (Mushotzky & Locwcnstein 1997; Boehringer et al. 1998; White 2000) 
and BeppoSAX (Ettori et al. 2001). Chandra showed that Abell 2390 is experiencing a minor 
merger event(see e.g. Allen et al. 2001; Vikhlinin et al. 2005; Baldi et al. 2007). Outside the 
cooling region, the average temperature is 11.5 keV (Allen et al. 2001). Abell 2390 has a 
small (less than 2' ) irregular radio halo, most likely related to the central AGN (hke the 
mini-halo of Perseus cluster, Bacchi et al. 2003). Its flux density is 63 mJ at 1.4 GHz and 
the equipartition magnetic fleld was estimated to be 1.3/iG (Bacchi et al. 2003). Using the 
model of the surface brightness profile of Abell 2390 determined with Chandra (e.g. see Allen 
et al. 2001, for details), we derive that selecting photons in XMM-Newton within a radius 
of 10' of the core includes virtually all cluster's emission. 



^The probability of observing a fluctuation has to be multiplied by the number of pixels in the BAT map 
(i.e. 2.9x10'^, Segreto et al. 2010). 
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Fig. 13. — BAT significance map of Abell 3667 witli superimposed X-ray contours from the 
ROSAT-PSPC (green) and radio 843 MHz SUMSS contours (white). 

There are two bright X-ray sources located within the selection region. Their coordi- 
nates are respectively R.A.(J2000)=21:53:40.7 and Decl.(J2000)=17:44:13.8 for the brightest 
source R.A. =21:53:34. 6 Decl. = 17:36:26. 8 for the dimmer one. Crawford et al. (2002) con- 
ducted follow-up observations of all sources detected with Chandra in the field of A2390. 
The sources reported above correspond to the sources A20 and A19 in their paper. Both 
these sources are AGN with A19 being a Seyfert 2 galaxy at z=0.305 and A19 a QSO at 
z= 1.6750. The spectrum of A20 is well fit by an absorbed power law with column den- 
sity of = 6.9+° X 10^1 cm-2 and a photon index of 1.52+°;^^ Its flux in the 2-10 keV 
band is 2.84+?|° X 10 erg cm s while the extrapolated flux to the 15-55 keV band is 
5.8x10"^^ erg cm^ s~^. The spectrum of A19 is compatible with an unabsorbed power-law 
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Fig. 14. — BAT spectrum of Abell 3667. The dashed hne is the best fit thermal model with 
a temperature of IS.T^g^/ keV. 

model with a photon index of 2.02l°:t^. Its 2-10 keV flux is 7.761^;^^ x 10"^^ erg cm-^ s"^ 
while extrapolated flux to the 15-55 keV band is 6.2x10"^^ erg cm^ s~^. It is clear that A19 
might contribute a non-negligible fraction (e.g. ~25 %) of the total flux detected in the BAT 
band while that is not the case for A20. Thus when analyzing the cluster emission (below) 
we will also include, in all spectral fits, an absorbed power-law component representing the 
spectrum A19. The parameters of this absorbed power-law will be allowed to vary within 
their 90 % CL reported above. 

A single-temperature plasma model (reported in Fig. 16) with a temperature of 9.471q 44 
and an abundance of 0.32±0.06 solar successfully fits the XMM-Newton and Swift/BAT 
data. Our results are in good agreement with those derived in the 0.5-40 keV band by 
BeppoSAX (Ettori et al. 2001). This fit is good (xVdof=409.9/375), but it leaves some 
residuals at high energy. We then tried adding a second thermal model and obtained a 
better flt (e.g. x^/dof=394. 5/372). The improvement in the is signiflcant and the F- 
test yields a probabihty of ~10~^ that it was produced by chance. The 'cold' and 'warm' 
components have a temperature of 3.76l^;g5keV and 13.08l2;69keV respectively, while their 
abundances are 0.461q;24 and 0.37lo;i4- Our results are in agreement with those obtained by 
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AUen et al. (2001) using Chandra. Indeed, they showed that the temperature of the plasma 
within 100 kpc of the core is <5keV while its temperature stays approximately constant 
at 11.5^16 keV beyond 200 kpc. Adding a power law model to the baseline thermal model 
improves the fit only marginally (Ax^=6.1 for 2 additional parameters), thus we consider 
the double-temperature thermal model as the best representation of our dataset. This fit, 
together with the single temperature thermal model, is shown in Fig. 16. 

The 99 % CL upper limit on the 50-100 keV non-thermal flux, estimated using a power- 
law with a photon index of 2.0, is 3.14x 10^^^ erg cm^^ s^^. In order to compute the lower 
limit on the intensity of the magnetic field (see Tab. 4), we adopt for this cluster a value of 
the spectral index of a=2. 

4. Clusters Magnetic Field 

The diffuse synchrotron radio emission (radio halos, relics and mini-halos) proves the 
existence of magnetic fields and relativistic electrons in the ICM. If the non-thermal X-ray 
emission results from IC scattering of the same population by the CMB, then the lack of a 
detection of a non-thermal component can be used to place a lower limit on the magnetic 
fields B in clusters. Indeed, the ratio of radio to IC flux scales proportionally to Fol- 
lowing Harris & Romanishin (1974), we estimate the lower limit on B (the volume averaged 
component along the line of sight) as explained in Ajello et al. (2009a), but taking into 
account the redshift correction. We model the IC emission as a power law with index 2 (see 
e.g. Reimer et al. 2004, for more details). The value of the diffuse radio flux is difficult to 
measure due to the presence of individual radio sources and to the variability of the spectral 
index with the distance from the center. Therefore the magnetic fleld intensities listed in 
Tab. 4 have to be taken as order of magnitude estimates. We flnd magnetic flelds that are 
typically a fraction of a /xG, thus far from equipartition. Note that A2390 is the only cluster 
for which we evaluate the magnetic field related to the radio mini-halo — and hence to the 
central AGN — rather than to a more extended radio halo or radio relic. It was not possible 
to estimate the lower limit on the magnetic field intensity for a few of the clusters reported in 
Tab. 4 for which there are no detections of radio-halos reported in the literature at this time. 
Given the fact that the sensitivity of BAT in its band is of the order of ~ 5 x 10^^^ erg cm~^ 
s~^ and thus comparable to the sensitivities reached (in other bands) by other observatories 
(e.g. HEAO-1, RXTE, Beppo-SAX etc.), the upper limits reported in Tab. 4 are similar to 
those obtained by other authors (e.g. see Rephaeli et al. 1987; Henriksen 1998; Rephaeli 
et al. 1999; Rossetti & Molendi 2004, and references therein). 
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5. Conclusions 

The present work combines Swift/BAT and XMM-Newton observations to investigate 
the presence of a hard X-ray excess in the spectra of 10 galaxy clusters detected in the 
ongoing BAT survey (Cusumano et al. 2010). Our results agree with our previous findings 
for a sample of 10 clusters (Ajello et al. 2009a) — i.e., most of the clusters' spectra are 
best described by a multi-temperature thermal model. The only exception is represented by 
the Bullet cluster and Abell 3667, for which we find evidence (at the 4.4(7 and 4.6(7 level 
respectively) for a hard X-ray excess. 

For the Bullet cluster, our data points to the existence of a power-law like component 
with a photon index of IMtlH and a 20-100 keV flux of 3.41^;^ x lO^^^erg cm-^ s^^ The 
flux of this component is found to be in good agreement with similar values reported by 
Petrosian et al. (2006) and Million & Allen (2009). Using the flux reported above and radio 
data available in the literature, we estimate that the volume average magnetic field should 
have an intensity of ~0.2 /iG. 

The case of Abell 3667 is different. Indeed for Abell 3667, the excess can be explained in 

terms of a hot component with a temperature of 13.5^2 2 keV. Our findings are in agreement 
with the results from Suzaku (Nakazawa et al. 2009). The lack of a central radio halo in 
Abell 3667, supports the thermal origin of the hard X-ray excess. 

The Norma cluster is a special cluster for a different reason. It is the second cluster, 
along with Coma (Ajello et al. 2009a), to be resolved spatially by BAT. The BAT spectrum 
shows (albeit with low statistics) that the temperature is around ~10keV and thus hotter 
than the temperature (~5keV) determined at lower energies with XMM-Newton. Since it 
is resolved by BAT a special care must be taken when analyzing data from this cluster and 
a detailed analysis will be presented elsewhere. 

Three of the detected clusters (PKS 0745-19, A1795 and A2390) have a bright cool core, 
while three (Bullet Cluster, A2256 and A3627) are undergoing a major merger. A1914 , A85 
and A3667 show signs of a minor merger as well. Six clusters (Abell 401, Bullet Cluster, 
PKS 0745-19, AbeU 1914, Abell 2256, and Abell 3667) have a radio-halo or a radio-relic. 
The best spectral fits are given by the sum of two thermal components for five clusters (A85, 
A401, PKS 0745-19, A3667 and A2390). The other four clusters (A1795, A1914, A2256, 
A2390) are successfully fit with a relatively hot single temperature profile. 

The upper hmit to the non-thermal emission (in the 50-100 keV band) is around 10~^^ erg 
cm~^^ s~^ for most clusters. Once again, our results indicate that the hard X-ray emission 
from galaxy clusters is mostly thermal and probably related to post-shock regions (in the 
case of merging clusters) or hot regions outside the cool core (in the case of relaxed clusters) . 
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It is reasonable to assume that the relativistic electrons observed in the radio band, produce 
a power-law like emission at higher energies due to IC scattering off CMB photons (e.g. 
Rephaeli 1979; Sarazin 1999). In this case, we used the presented results to obtain lower 
limits on the clusters' magnetic fields. In this way, we find magnetic fields of the order of a 
fraction of yuGauss, that are far from equipartition and in agreement with previous similar 
estimates. These limits are generally a factor 10 below the estimates obtained using Faraday 
rotation measures (e.g. Clarke et al. 2001; Guidetti et al. 2008; Bonafede et al. 2010). 

The Bullet cluster is the only one that stands out among the clusters detected so far 
by BAT for the evidence of a power-law like, hard X-ray excess. However, many factors 
make this cluster special and unique; among them: the intermediate redshift and the violent 
merging activity. The merging process powers shocks (e.g. Markevitch et al. 2002) where 
CRs can be accelerated efficiently. The energy density of the CMB, whose photons constitute 
the targets for the electrons, scales with (1 -|- z^, and thus is a factor ~2.8 larger than at 
redshift zero. Both things probably concur in producing the 'bright' non-thermal component 
observed in this cluster. 
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Fig. 15. — XMM-Newton and BAT data for Abell 3667 fitted with: 1) a single thermal 
model (top), 2) the sum of two thermal models (middle), and 3) the sum of a thermal and 
a power-law model. 
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Fig. 16. — Left Panel: Spectrum of Abell 2390 fitted with a single thermal model plus an 

absorbed power law (dashed line) to account for the emission of the AGN A19 (see text for 
details). Right Panel: Spectrum of Abell 2390 fitted with the sum of two thermal models 
(thin solid and short dashed line) and an absorbed power law (dashed line) to account for 
the emission of the AGN A19 (see text for details). 



Table 3. Spectral Properties of BAT detected Galaxy Clusters (errors are 90% C.L.) 



NAME 


z 


Flux^ 
(10-12 j,gs) 


(10"*^ erg s-i) 


kT 

(keV) 


r/kT 


model 


xVdof 


Abell 85 


0.0521 


r ic+O.Sl 


q Ol +0.52 

'J-°l-0.82 


6.09l°:g 


1 79+0.32 
1- ' ^-0.06 


apec + apec 


602.1/619 


Abell 401 


0.074 


6.39^°|i 




8.6iin^ 


Z.UO_Q 45 


apec + apec 


732.4/652 


Bullet 


0.296 


5.iol?;fo 


176+^1 


14.77tJ:^^ 


1-861S:?! 


apec + pow 


501.7/511 


PKS 0745-19 


0.103 


6.93tof9 


23.2111 


7.96^°:g 




apec + apec 


587.1/578 


Abell 1795 


0.062 


2 n5+°-i^ 


r) o'7~l~0.23 
-0.20 


4.82±°:i; 




apec 


892.1/1275 


Abell 1914 


0.171 


4.29tl;°^ 


46.6+:^-^ 






apec 


355.1/351 


Abell 2256 


0.0581 




4.04tl:i? 


8.84t°J? 




apec 


445.6/445 


Abell 3627'' 


0.0168 










brcm 


14.7/14 


Abell 3667 


0.0556 




5.65t}°^ 




13.5^^:1 


apec + apec 


569.8/542 


Abell 2390 


0.231 


9 1 q+0.26 
^•-'^''-0.25 


52.5tl^ 


13.08l2;69 


q 7^+2.80 
"^•"^-1.61 


apec + apec 


394.5/372 



*Flux and Luminosities are computed in the 15-55 keV band. 
^For this cluster only BAT data were used. 
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Table 4. Non-thermal emission from combined XMM-Newton and BAT data. 



NAME 


i^50-100keV 


B*^ 




(10-12 erg cm2 s"!) 


(mG) 


Abell 85 


<2.51 


~ 0.6 


Abell 401 


<0.22 


~ 0.4 


Bullet 


1 ro+0.43 
-l-.Of5_o.47 


~ 0.16 


PKS 0745-19 


<1.6 


~0.5 


Abell 1795 


<1.38 


/ 


Abell 1914 


<1.08 


~ 0.3 


Abell 2256 


<0.19 


~ 0.6 


Abell 3667 


2 98+'^-^'^ 

Z.»0_o.73 


/ 


Abell 2390 


<0.25 


~ 0.8 



^The flux has been estimated using a 
power-law spectrum with a photon index 
of 2.0 in the 1-200 keV energy band. Up- 
per limits are 99% CL while errors are 
90% CL. 

''In order to compute the intensity of the 
magnetic field we used the radio data listed 
in Section 2. When a was not available, we 
adopted a — 2. 



